MIT-CTP-3404 



m 
o 
o 

(N 



The self-consistent bounce: an improved nucleation rate 

Yoav BergneiQ 

Center for Theoretical Physics, Massachusetts Institute of Technology, 
77 Massachusetts Avenue, Cambridge MA 02139 

Luis M. A. BettencourtQ 

Los Alamos National Laboratory, MS B256, Los Alamos NM 87545 



O ; (Dated: February 1, 2008) 

! Abstract 

CN ' We generalize the standard computation of homogeneous nucleation theory at zero temperature 

. to a scenario in which the bubble shape is determined self-consistently with its quantum fluc- 

', tuations. Studying two scalar models in 1+1 dimensions, we find the self-consistent bounce by 

I employing a two-particle irreducible (2PI) effective action in imaginary time at the level of the 

o ' 

rate exponent. We use collective coordinates to account for the translational invariance and the 

Q , growth instability of the bubble and finally present a new nucleation rate prefactor. We com- 



Hartree approximation. We thus obtain an effective single bounce action which determines the 



pare the results with those obtained using the standard 1-loop approximation and show that the 



^ ■ self-consistent rate can differ by several orders of magnitude. 
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I. INTRODUCTION 



During first order phase transitions, the macroscopic state of a system changes suddenly 
from the metastable phase (or false vacuum) to the stable phase (or true vacuum). If the 
transition is strongly first order it is expected to proceed by the spontaneous creation of 
spherically symmetric droplets - or bubbles - of the stable phase which, if large enough, can 
grow to consume the metastable phase. 

The modern theory of bubble nucleation is due principally to Langer who was con- 
cerned with models of statistical physics and to Voloshin ^ and Coleman and Callan P, 0] 
who developed similar ideas in the context of relativistic (zero temperature) quantum field 
theory. The unification of Langer and Coleman's approach to finite temperature quantum 
field theory was later obtained by Affleck and Linde 0], among others. 

Although the departure point for the theory of nucleation is completely general, in prac- 
tice the nucleation rate is computed almost without exception in the semiclassical approxi- 
mation, including fluctuations only at most to 1-loop order. Even for semirealistic models, 
these lowest order calculations can be computationally nontrivial . We quote the now well- 
known zero temperature result of Ref. 0| that for D spacetime dimensions, the nucleation 
rate per unit volume T/V is given by 



T/V = Ae-^, B = S,i[ct>,]/h, A={ — 



det' [-d + V"{<Pb)\ 



det [-d + V"{(j)+)] 



'1/2 

(1) 



where Sd = J d^x{dx4>Y /2 + V{(j)) is the classical Euchdean action of the "bounce" and det' 
denotes, as usual, the fluctuation determinant with D zero modes excluded. 

Nevertheless in certain instances, an improvement on the 1-loop computation is neces- 
sary, for example when the flrst order character of the transition is itself due to radiative 
corrections or to the presence of other flelds, as in Higgs + gauge theories. It would seem 
that, at least to lowest order again, these effects might be accounted for by substituting an 
effective potential for the classical one in the bounce computation or more generally em- 
ploying an effective action . But this is not a simple task as fluctuations exist over the 
bubble background, and this background may distort significantly the properties of the low 
lying spectrum relative to its perturbative form. 

One resolution to this problem is a different treatment of low and high energy fluctua- 



tions, for example in the context of a coarse-grained effective action j9|, |lO|, lUl- In such a 



nign 

m 



procedure, however, much care must be taken to avoid double counting the effects of fluctu- 
ations, i.e. in their accounting for distortion of the bubble while also providing corrections 
to the background in which the bubble exists. Certain physical effects such as fluctuation 
backreaction on the classical bubble profile are usually absent from these calculations. Fi- 
nally, lattice Monte Carlo methods have also been used to st udy nucleation rates, including 
more recently for the case of radiatively induced transitions |l2l Il3l|. 

In this paper we propose a self-consistently improved nucleation rate by including all of 
the above features in a semi-analytical computation We shall study two scalar models in 
1+1 dimensions and include the self-coupling of fluctuations as well as their backreaction on 
the mean field nonperturbatively through the use of a two-particle irreducible (2PI) effective 
action formalism. We have also used this technique recently to calculate the self-consistent 
quantum energies of topological defects 

The organization of the paper is as follows: in Sec.|n]we review aspects of the constuction 
of the rate as defined in Eq. (0) as some of these will need to be modified later. We define 
our scalar models in Sec. IIIII and give a short overview of the 2PI effective action formalism 
we employ. We also define the quantum bounce in terms of coupled Euclidean spacetime 
equations for the mean field and the coincident two-point function. In Sec. lIVI we discuss the 
renormalization structure of counterterms necessary to make these calculations finite. We 
present numerical results for the quantum bounce profiles and corresponding effective action 
calculations in SecEJand also discuss the spectrum of fluctuations in the self-consistent case 
vs at 1-loop order. Finally in Sec. I VII we discuss how to extract the nucleation probability in 
the self-consistent case and propose a new nucleation rate. A brief discussion of the results 
is offered in Sec. IVIII 



II. STRUCTURE OF THE RATE CALCULATION 

The 1-loop approximation to the rate calculation leading to Eq. has several special 
features, which we now highlight in order to be able to contrast them with our results below. 
We will take h = 1 throughout and work in /^-dimensional Euclidean spacetime, although 
later we shall specialize to D = 2. 



The starting point of any rate calculation is the transition amplitude 

{M^)\Mx'))^M J ci[0]e-^[*l/^ = Z, (2) 

with appropriate boundary conditions on the fields. We are interested in a transition am- 
plitude between the false and true vacua, i.e. the vacuum expectation value of the field is 
given by the true vacuum at some spacetime point, taken to be the origin {(^{x'^ — 0)) — (p^, 
and in the false vacuum at spacetime infinity {(f){x^ — > oo)) = 0+. The normalization J\f 
ensures that 

^ {cP4xM4x')) ^ 1. (3) 

In the semiclassical approximation the path integral is evaluated by expanding around the 
field configurations that extremize the classical action, subject to these boundary conditions. 
The simplest such configuration is the classical bounce Well separated multi-bounce 
configurations also qualify and need to be summed over among the action extrema. 

The single bounce path integral can be evaluated semiclassically (i.e. to 1-loop order) by 
expanding the quantum field in linearized fluctuations around the classical bounce configu- 
ration 

= 06 + X] CtnV'n- (4) 
n 

Then the (renormalized) path integral becomes 



Z = J\fAe-^; 



det [-d + V"i(l)b)] 



-1/2 



B = S^[<j>t,] (5) 



det [-d+V"{(j)+] 

where we took the fluctuations to be eigenstates of the linearized operator G^^ — —d+V" {(/)). 
Deflned as such, the fluctuations are non self-interacting and do not backreact on the proflle 
0b (i.e. neither the operator Go nor the equation of motion solved by 0^ depends on the 

fluctuations ipn)- 

Note that Z is dimensionless as it should be for a probability. There is a simple functional 
relation between Z and the effective action Feff. In the absence of external sources 

Z = e-^«ff (6) 

and the 1-loop effective action takes the well known form Feff — B — ^InA. 

Given the transition amplitude due to one bounce, it is straightforward to generalize to 
any number of bounces. This is the sum of the 1-bounce contribution plus the 2-bounce 



contribution, which for well separated configurations becomes the product of two single 
bounces, etc. The result of resumming this series is the exponentiation of the single bounce 
amplitude i.e. 

Z = exp[e-^A] = expfe-'^-^s] (7) 

The nucleation rate (per unit volume) has dimensions of inverse spacetime volume and 
can be obtained directly from the multi-bounce amplitude as 

T/V = 2lm(^^\nZ^ , (8) 

where fl is the spacetime volume. (We regret the confusing notation which arises since the 
capital letter F serves as the conventional standard for both the rate and the effective action. 
We use Fefi where necessary to avoid ambiguity.) 

To see that this result is finite in the infinite volume limit, we must analyze the deter- 
minant ratio in the prefactor A in greater detail. To 1-loop order, it is well known that 
the fiuctuation spectrum has one negative eigenvalue (an imaginary frequency) and D zero 
eigenvalues. The former keeps track of the instability of the false vacuum. In a theory 
without backreaction this instability is connected with a nonconservation of a probability 
current j^lsj, a mode that will grow unchecked in the large spacetime volume. The D zero 
modes are proportional to the D variations dg(j)b, and represent the translational invariance 
of the "center" of the classical bounce. 

Each eigenvalue in the determinant A results from a Gaussian integration 

l^2-l,,e-*i*. (9) 

This integral is clearly divergent if u"^ is negative, but it is rendered sensible by deformation 
of the contour into the complex plane (in effect taking uj"^ I^^D- Thus the amplitude 
Z is imaginary, reflecting the fluctuation instability. (As the nucleation-or decay-rate is 
determined by the imaginary part of the transition amplitude, one simply excises the factor 
of i). 

In the case a; — the corresponding integral is of course not Gaussian at all. It becomes, 
up to a Jacobian factor, the volume and diverges with it. This is clear from the change of 
variables 

1/2 



[d(f)\i = Jdz, J 



d^xi^^ 
ox 



(10) 



lim — = \ — J 

^-^0 uj V 27r 



dz 



J 



27T 



(11) 



with the hnear volume n^^^ oo. What we have done is exchange the functional integration 
over the field direction [d(f)] for an integration over a collective coordinate z. The zero modes 
thus give rise to an infrared divergence (related to the limit of infinite spacetime volume) 
in the determinant pref actor. This factor of volume is happily canceled as the probability 
amplitude, taken per unit spacetime volume, becomes the desired rate per volume of Eq. ((H). 
Namely 



1 



det [-d + U"i^b)] 



det [-d + U"{(j)+)] 



-1/2 



J 



D 



det' [~d + U"i(Ph)] 



det [-d+U"{(j)+)] 



-1/2 



(12) 



(In a finite volume application, it is quite reasonable that the nucleation rate grows as a 
function of volume.) 

There are two distinct ways in which this standard calculation is approximate. The 
first is the computation of the single-bounce effective action. The semi-classical or 1-loop 
approximation is a lowest order calculation of the effects of fluctuations. It completely 
neglects both fluctuation self-interactions and their backreaction on the classical fleld proflle. 
Because barrier nucleation is a Boltzmann-suppressed process, exponentially sensitive to 
corrections, the inclusion of these effects may lead to large quantitative differences. A 
computation of the nucleation rate that includes these effects is the subject of the present 
paper. 

The second approximation involved in the rate calculation is the "dilute instanton ap- 
proximation" J, 31 1 for the multi-bounce conflgurations. This is in principle justifled so long 
as the nucleation rate per volume is sufficiently small. It will not be a good approximation 
though in cases where nucleation may be enhanced by earlier bubbles or by the presence 
of other nonperturbative fluctuations over the homogeneous background. Although some 
of these situations could be addressed with the techniques developed here we will not ex- 
plore them in the present manuscript. We will thus continue to assume the dilute instanton 
approximation as an unaltered ingredient of the rate computation. 



III. 2PI EFFECTIVE ACTION 



The traditional starting point for the calculation of the single bounce is the classical 

n 

action. We refer the reader to the Ref. y| for that formalism and do not summarize it 



extensively here. Suffice it to say that the classical bounce is an 0(i5)-symmetric solution 
to the Euclidean field equations in D spacetime dimensions satisfying certain boundary 
conditions, in particular that the field assumes a true vacuum expectation value in a local 
region around the origin and a false vacuum expectation value far away and to infinity. In our 
development of the self-consistent bounce, we begin instead with the two-particle irreducible 
(2PI) effective action for the field and the two-point function G{x, x') developed in the 
relativistic context by Cornwall, Jackiw and Tomboulis (CJT) jlS^ . 
The effective action functional is given by 

r(0, G) = + ^zTr InZ^oG-' + ^zTrD-^(0)G + 12(0, G) + const (13) 

In this notation we have the classical action 

/(0) = I d'xc[<p], m = ld,<pdy - v(0) (14) 

and the operator 

Dq is the free propagator. Finally r2 sums the 2PI vacuum to vacuum diagrams which 
furnish an expansion in the number of loops or, equivalently, in powers of h [h = c = 1 below 
and throughout). Stationarizing the action leads to the equations: 

61 1 6D-^ 5V2 , 
We shall consider two symmetry-breaking forms for the scalar potential 



V^6(0) = ^A02 



Xv 



a 



= ^ (18) 



Xa 

These potentials represent more or less generic models with true and false vacua at the 
classical level and have been cast in the form that the potential difference between the true 
and false vacua is e. The symmetry is explicitly broken in the quartic model while a — — 
discrete symmetry still exists in the sextic model (where the false vacuum is the symmetric 



vertex 



cubic -Xcj) -2A(503 - 2a^(l)) 
quartic -| -|(1502 _ 2^2) 

quintic — 3A0 

sextic — ^ 

TABLE I: Interaction vertices used in the construction of the 2PI effective action for two models 
minimum = 0). The quantization of solitons in ( 1 + 1 ) -dimensional models of precisely 



these types with e = were originally examined in Ref. 



1( 



For these models we obtain the following equations: for V4, 



□0+-A 



5To 

+ 3A0G(x,x) - -f = 



iG-\x, x') + (n^. - + 3A02) 5^(x -x')-2 



5V2 
5G 



(20) 
(21) 



and for Ve, 

□0 + A0 



2e 
A^ 



6A (2aV-50^)G(x,x) 







^(x,x') + 



+ A(a^ 



2e 
A^ 



12Aa>^ + 15A0' 



6'^{x-x')-2 



5T2 
5G 



(22) 
(23) 



Up to this point we have made no approximations, but we have left r2 as the sum of 
2PI vacuum-to-vacuum diagrams with interaction vertices determined by the shifted La- 
grangian and lines representing the full propagator G. For our models, we can summarize 
the interactions in Table HI 

It is of practical necessity in the computation to select only certain 2PI diagrams, and such 
a choice will correspond to a truncation in the Dyson-Schwinger hierarchy of equations for 
the correlation functions of the model. It is straightforward to write down the full evolution 
equations in the Hartree approximation, which amounts to including only the 2PI "bubble" 
diagrams shown in Fig. ^ In 0{N) models, this truncation is similar to the leading-order 
large- approximation. While the Hartree approximation is not a controlled, systematic ex- 
pansion, it is understood to be equivalent to a Gaussian variational ansatz in the Schrodinger 
functional formalism [2^ and results in Hamiltonian dynamics j2^. Furthermore, a study 



FIG. 1: Vacuum bubble diagrams in the Hartree approximation for the 2PI effective action of 
our models. For the classical action V4, only the left diagram contributes. Lines denote the full 
propagator G in the nontrivial <j) background. 



of the quantum energy of solitons showed good agreement between the Hartree values and 



"exact" lattice Monte Carlo results 



16 



T^^ = -^A j d^xG{x,xf 



17|. Proceeding with the approximation, we find 

(24) 



J- 2 



■^A f d^x{15(l){x)'' -2a^)G{x,xf -^X / d^xG{x,x)^ 



so we obtain for V4 (after inverting the G equation) 



□0+-A 



2(^((f)^-v^)-— +3X(f)G(x,x) =0 



[n^ - Xv^ + 3A02 + 3XG{x, x)] G{x, x') = -t5^{x - x') 



and for Vk. 



□ + A 



6A {2a^ - 502) G(x, x) + 45AG(x, xf\ = 



□x + A(a^ 



2e 
A^ 



12Aa^0^ + 15A0^ 



+6A(1502 - 2c?)G{x, x) + 45AG(x, xf] G{x, x') = -i5\x - x') 



(25) 

(26) 
(27) 



(28) 



(29) 



This formal procedure has absorbed the approximate quantum dynamics of the scalar field 
into a truncated set of equations to be solved simultaneously. The two-point function 
G{x, x') now appears explicitly as a "sort of" Green's function, albeit of a nonlinear operator 
which depends on its value at coincident spacetime points. Typically we would proceed by 
decomposing G{x, x') in a mode basis (related to the decomposition of the field operator of 
which it is a correlation function). The mode functions and the field then satisfy a tower 



of partial differential equations which can be stepped forward in time on a computer. This 
would be a prescription for the nonequilibrium evolution given Cauchy-type initial conditions 
221. 



A. Imaginary Time 

We shall alternately consider attempting to solve Eqs. ()26|) - ()29p in imaginary time. Pass- 
ing to Euclidean spacetime and framing the partial differential equations as elliptical equa- 
tions instead of hyperbolic ones does not in and of itself simplify the problem unless we can 
exploit special symmetries. This is precisely what is done in the case of the bounce, where 
as a consequence of Lorentz invariance we seek a rotationally invariant classical solution, 
i.e. we need solve only for radial functions. We can now see why at least at the level of the 
Hartree approximation, we can attempt a similar sleight of hand on the two-point function 
and the field simultaneously. Wherever the two-point function appears in the evolution op- 
erators, it appears with coincident spacetime arguments. Therefore, we can seek solutions to 
Eqs. ()26|l - (j29|l in imaginary time in which both and G{x,x) depend only on the radial 
coordinate in Euclidean spacetime. 

To wit, the imaginary time formulation of the problem entails, 

t^iT, o^Oe = -W^ (30) 

or going to polar coordinates 

1 d f d\ 1 

"""-—pd-pypd-pj-y^w ^^^^ 

We seek a solution to the field VEV and two-point function equations where the field VEV 
and the coincident two-point function are purely radial, i.e. 

= 0(p), G{p,e;p,9) = G{p) 

while in general there remains ^-dependence in G{p,6; p',9'). It is consistent with these 
constraints to separate variables in which case G takes the form 

G(p, 9- p', 9') = aniRni{p)Rni{py^'~''^ (32) 

n,l 



The (f) equations and pKjl become 



-A 



V- - 



+ 3A0G(p) = 
2e 



— a 



6A (2a2 - 502) + Ah\G{pf] <i){p) = 



and the equations 



[□e + x{p)] G{p, e- p', e') = -6{p - p')6{d - e' 

X4(p) = -m2 + 3A02 + 3AG(p). 



(33) 



(34) 



(35) 



2e 
A^ 



12Aa>' + 15A0^ + 6A(150' - 2a')G(p) + 45AG(p)'. 



lead to the following "nonlinear eigenvalue" problem for the radial modes 



1 d 



d 



pdp\dp)^p^^ ^^^^ 



Rr, 



with the identification 



''nl 



nil 



Rliip) 



n,l 



provided the orthonormality and integrability condition holds for the Rni{p), 



2-K / pdpRnl{p)Rmk{p) = S 



nm,kl 



(36) 



(37) 



(38) 



We must also impose boundary conditions at the origin p = 0, where for analyticity 
we require that G have zero derivative, and at large values of p where G should approach 
its value in the false vacuum. We impose the first condition on the radial mode functions 
Rni{p) whereby for / = the functions must have zero derivative at the origin and for / 7^ 
they must vanish. The condition at large p can be made consistent with Dirichlet boundary 
conditions. This latter condition is especially suitable if the renormalized value of G is 
fixed to be zero in the false vacuum. We shall discuss this renormalization condition in the 
following section. 



IV. RENORMALIZATION OF SELF-CONSISTENT EFFECTIVE ACTION 



Although any scalar model in 1+1 dimensions is super-renormalizable, details of the 
renormalization differ slightly in the self-consistent effective action from the usual loop ex- 
pansion. Moreover numerical implementation of renormalization schemes comes with its 



own subtleties. Here we detail the procedure for a simplified 0^ model, i.e. with no sym- 
metry breaking parameter e, in two steps. We first demonstrate the analytical cancellation 
of divergences assuming a homogeneous background. This may be considered a review of 
the renormalization of the self-consistent effective potential as worked out previously, for 
example in Ref . |23j . Next we describe the changes which occur for the inhomogeneous case 
of interest and for a finite volume with a discrete spectrum-as is the case in any numerical 
analysis. 

For homogeneous fields we consider the effective potential and pass to Euclidean momen- 
tum space for the translationally invariant two-point function, i.e. 

V"(0,G) [ (fx = T{^,G), (39) 



such that 



2 



V{<P,G) = V,,{^) + ^ IjnDoG-' + ^ Ijd-'G -1] + ^[I^G{k)\~ (41) 
Now setting |^ = 0, we have 

G-\k) = D-\k) + 3\ I G{k) (42) 

Jk 

and taking as an ansatz G{k) = pqr^, we obtain a gap equation for M^: 

M= = -»=+3A^= + 3A/|ijj^]^ (43) 

The integral in the gap equation is logarithmically divergent, hence is so far undefined. 
We now demonstrate the counterterm procedure which renormalizes this equation and also 
renormalizes the effective potential. We first rewrite the equation as 

Af = -,n^ + 3A^^ + 3A (/ (g,^,^]^ - /) + 3A/ (44) 

and impose a momentum cutoff A with 

^ n t d'k 1 1 A2 



(27r)2 ' J (27r)2 F + /ig 47r fi^ 



This implements the regularization of the integrals in Eq. ()44j) . The new scale /io is necessary 
to prevent an infrared logarithmic divergence in the infinite volume limit. The renormaliza- 
tion procedure is now carried out by identifying the bare parameter as a function of the 
cutoffs and defining a renormalized (physical) mass by 

-m% = -m^ + 5m\ = -m^ + 3A/ (46) 

In terms of the renormalized mass, the gap equation is finite 

M^ = -m| + 3A0^ + |ln^ (47) 
= -m^ + 3A02 + 3AG/(M). (48) 

and we can see that by a choice of renormalization scale (/ip = 2m|.), the self-consistent 
fluctuation mass in the classical vacuum 0q = fn^/X can be fixed to its "classical" value 
= 2m|.. We shall consider fiQ to be a physical scale, a function of the physical mass and 
couplings. 

We now express the Hartree resummed effective potential in terms of its contributions at 
different orders, following Ref. Q]. 

Veff = Vo + Vi + V2 (49) 

Substituting Eq. (j42|) into (jUj), we now have 

1 



V, = iTrlnDoG-^ = ^ j 

















(2vr) 


3A 




~T 


(/ 



2 

1 r d:^k 

In 



2 



3A^, ,2 3X f f cPk 1 V , . 

Note that as a simplification this classical potential has no symmetry breaking term and no 
overall constant. The constant has no physical significance, but we will in any case make 
the shift in the potential explicit below by subtracting the value at a specified vacuum. The 
symmetry breaking coupling introduces no new infinities and requires no renormalization. 
We use (HSl) to write 

G(x, x) = ^ (M^ - - 3A02) , (51) 
3A 



and 

Vo + V, = -^- 1a0^ + ^M'm' + lM'<f>\ (52) 

(where an additive constant has again been removed). We rewrite in terms of renormahzed 
parameters, 



1 . A' 



—m 



rrij^ — SXI, I = Gq{x,x) = —In— (53) 



471 

and use the renormahzed gap equation (jlSj) to write 



such that 



Im^G^M) = + ^ - lM20^ (54) 



1 1 1 

iVo + V,)^ = — - -A0^ - -M'G,{M) - -MH. (55) 



Finahy, the one-loop contribution to the potential comes to 



(56) 



2 J (27r)2 + 



After shifting the combined tree-level and 2PI pieces of the effective potential by their 
contribution in the classical vacuum (p^ = mj^/X (where the renormalization scale is set), 
the renormahzed total is 

Veff = -IX{<P' - 4) + + liM' - f^ll (58) 

which is finite. 

The generalization of the effective potential to the effective Euclidean action for inho- 
mogeneous mean field background now follows. We no longer have a simple gap equation, 
however the classical field x defined in Eq. fl35|) now replaces the in the above analysis. 
The counterterm renormalization goes through and renders the tree level and 2PI compo- 
nents as a finite piece, expressed in terms of renormahzed parameters, plus a counterterm: 



0+2 



ro+2 -^I j dx (x(x) - xo) (59) 



where 



= I d'x( i(V0)^ + ^ - - lx{x)G,{x) - ^ + ^-X4>t ] . (60) 



We have again chosen the mass scale of the counterterm to set G/(0o) = 0, so 

X{x) = -ml + 3X(f)'^{x) + 3XGf{x), 
Xo = -mji + 3X(f)l, (61) 

and we could have equivalently written 

= / d'x - K/(0o) - ^Gfix,xf^ , (62) 



with Vci expressed in terms of renormalized (finite) parameters. Finally, the one-loop order 
contribution in the inhomogeneous case is obtained from the eigenvalues of Eq. ()36p . 

by which we mean that are the eigenvalues in the presence of the spatially varying 
field x(p) and (c<Jq)„/ are the eigenvalues in the vacuum defined by Xo- This term is still 
logarithmically divergent, but this divergence is cancelled by the second term in Eq. ()59|l 
with 

The renormalization of the 0^ model follows along the same lines as above, although it 
is somewhat more cumbersome. In that case there occurs a renormalization of both the 
quadratic and quartic couplings. Both divergences have the same origin in the logarithmic 
loop integral of a gap equation analogous to Eq. PHj) . i.e. there is still only one renormal- 
ization scale. 

There remains a subtlety in the practical computation of this one-loop term (jHS)), which 
has been studied in the context of quantum corrections to static solitons j^]. In that case, 
the quantity of interest is the sum over zero-point energies, or more precisely the difference 
of this sum in the presence of a soliton from this sum in the vacuum. Rebhan and van 
Nieuwenhuizen have pointed out that in the lattice discretized sum (for example in a finite 
volume) care must be taken in implementing a frequency cutoff for the sum such that the 
same number of modes are included in the vacuum sum and in the nontrivial sum (ambiguity 
arises because of both the presence of zero modes at the one-loop level and the phase shifts 
of the "continuum" states). We shall refer to this as a frequency cutoff with mode parity 



(FCMP) instead of a mode number cutoff as in Ref. 



. 3- 



The analysis here of a spherically 



symmetric background in higher dimensions must follow the same prescription, and moreover 
the FCMP must be implemented separately for each angular quantum number. In other 
words, there is one frequency cutoff defined by the lattice spacing, but this leads to the 
inclusion of a different number of radial modes of each angular index included in the sum. 



V. THE SINGLE BOUNCE EFFECTIVE ACTION 

We have implemented the procedure described in the previous sections on a desktop 
workstation. Namely, we compute the vacuum counterterm (jMj) by solving the radial eigen- 
value problem (j36|l in the false vacuum of each model. We then introduce a nontrivial mean 
field bubble of true vacuum (using the semi-classical bounce as an initial guess^) and solve 
the renormalized versions of Eqs. simultaneously by successive iteration. We use a 

standard relaxation routine 25| to solve Eq. ()33|). together with a set of LAPACK packages 
3| to solve the eigenvalue problem Eq. ljH^ . At each step, the fluctuation eigenproblem 
is solved to produce a new intermediate Gi^ti^), which we combine with the old value to 
produce a new trial Gncw = (1 — 7)G'oid + jGint, until convergence in ip and G{x) is reached, 
up to a specified precision. The adjustable parameter < 7 < 1 controls the size of the 
update and can be adapted to optimize convergence. 

This numerical procedure will find the correct shape of the mean-field profile and the 
coincident two-point function if the radius is constrained. Unconstrained, however, it may 
fail to find the correct radius of the bubble, given the initial guess, since the stationary point 
of the action is a saddle point, not a minimum. (In other words, the successive iteration will 
not flow towards the saddle point.) The self-consistent bounce typically has a smaller radius 
than the classical bounce. We effectively constrain the radius first at its classical value and 
then vary the radius to find the actual saddle point. 

The action is computed at each iterative step. It should be noted that the very first step 
of the procedure amounts to computing the fluctuation contribution to the action in the 
presence of the classical bounce without backreaction on the mean field or self-coupling of 



^ The initial guess is highly arbitrary and our use of the semi-classical bounce does not for example preclude 
the application of this method to a theory with a radiatively induced phase transition. The nature of 
the algorithm allows our initial guess to relax into the correct bounce profile provided the boundary 
conditions, e.g. the approximate symmetry breaking vacuum expectation values, are initially satisfied. 
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FIG. 2: The self-consistent bounce and coincident two-point function for the cf)^ model with 
parameters \/ ^j? = 0.03, e//x^ = 0.15 are shown in solid lines. The dashed line is the (semi) classical 
bounce for this model. (Note that the y-axis for G{p) is shifted with respect to the axis labeled by 
0(/9); G still goes to zero in the false vacuum.) Inset: Euclidean action for various approximations. 

the fluctuation modes. In other words it is the standard 1-loop correction. Thus we are 
able to compare the value of the Euclidean action at the classical, 1-loop and self-consistent 
levels^ . 

In Figs. I21IS1 we show the classical bounce along with the self-consistent bounce and the 
coincident two-point function for our two models. In the inset of each we give the value 
of the Euclidean action at the different orders of approximation. As the figures indicate, 
it is possible for the self-consistent bounce profile to differ significantly, both in radius and 
thickness, from the semi-classical prediction. Moreover the action, which will enter into our 
rate calculation, is dramatically different. 

It is important to mention also the difference between the self-consistent fluctuation 



^ In the one loop case we excise the zero modes and take the absolute value of the negative mode, as usual. 
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FIG. 3: The self-consistent bounce and coincident two-point function for the (/>^ model with 
parameters X/ ^j? = l,e//i^ = 0.12 are shown in solid lines. The dashed line is the (semi) classical 
bounce for this model. Inset: Euclidean action for various approximations. 

spectrum and the 1-loop fluctuation spectrum. We display a plot of the low-lying spectra 
for the (f)^ model in Fig. EJ Modes are labeled by radial and angular quantum numbers n and 
/ in Eq.()36p. and in the figure the x-axis proceeds along the increments in the angular number 
/. The 1-loop and self-consistent spectra for each / are shown side-by-side for comparison. 
The numerical spectrum computed in a finite volume is of course discrete, but the lowest 
four states shown are actually bound states. 

What is noteworthy is that the negative mode and the two zero modes of the 1-loop 
spectrum (each mode for Z 7^ is doubly degenerate) are no longer negative or zero in the 
self-consistent spectrum. In fact, this is hardly surprising: once the fluctuation modes are 
allowed to interact in the presence of the bounce background, the spectrum will of course 
be shifted to positive values as the modes stabilize each other jisl . |Q. Needless to say 
this does not mean that translational invariance of the center of the bounce has been lost. 
Translational invariance is however no longer manifest at the level of linearized fluctuations 
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FIG. 4: Spectra of the low lying fluctuation modes for the (j) model. Plotted are the first few 
modes in the discrete spectrum for the first few partial waves, labeled by the quantum number 
I. The 1-loop spectra are shown side-by-side with the self-consistent spectra, with labels on the 
X-axis indicating IL for 1-loop and SC for self-consistent. 

since the mean field now cannot be translated without also translating the coincident two- 
point function G. We will return to the issues raised by this fact in the next section. For 
now it is perhaps most clearly illustrated at the level of the equations of motion. 

Consider for example, the 0^ model. At the 1-loop level, we start with the usual equation 
of motion for the field and for linearized fluctuations around the bounce. 

□0 + 1/'(0) = O (65) 
[a + v"i^)]M^)=^lM^) (66) 

It is evident by taking the derivative of the field equation that each of the modes correspond- 
ing to translation of the field satisfy the fluctuation equation with zero eigenvalue. A 
similar trick will not work on the set of self-consistent equations in ()26|) and ()27p . which we 



recast below as follows: 



+ V'{(j)) + 3X(f)G{x, x) = (67) 
[□ + + 3AG(x, x)] M^) = ^IM^) (68) 

The coincident two-point function G{x, x) now acts as an effective source in the fluctuation 
equation. 

VI. COLLECTIVE COORDINATES AND THE NUCLEATION RATE 

As we have just seen the low lying fluctuation spectrum in the self-consistent case is 
qualitatively different from the spectrum at 1-loop. The absence of the negative and zero 
modes is a direct consequence of the mode self-consistency and now implies that translational 
invariance is no longer a symmetry of the mean-field bounce profile. Nevertheless it clearly 
remains a symmetry of the transition amplitude. This situation is reminiscent of 1-loop 
calculations where the path inte gra l possesses more (e.g. internal) symmetry than exhibited 



by the bounce solution j27|, |28|, |29|. In these cases the contributions from such degrees of 
freedom must be accounted for in the transition amplitude (and consequently the rate of 
decay) via collective coordinates. That is the goal of the present section. 

The existence of one negative eigenmode and D zero eigenmodes in the 1-loop case carries 
important physical information about the nature of the bounce solution and its symmetries. 
The former indicates that there is an instability related with the bubble's expansion or 
contraction, whereas the latter are a consequence of the translational invariance of the bubble 
as a whole. In 1-1-1 dimensions these modes are related as they refer to displacements of the 
bubble walls in opposite directions or in the same direction, respectively. The translational 
and dilatational excitations can in fact be built up from the zero modes of the bubble walls 
as their symmetric and their anti-symmetric combinations in the thin wall approximation 

In the self-consistent case the bubble growth instability and its translational invari- 
ance persist globally in the effective action, but are no longer manifest at the level of 
the fluctuation spectrum, which now consists solely of positive, "vibrational" modes. As 
we have alluded, translational invariance now requires simultaneous displacements of both 
(p{x) — > (p{x + a) and G{x, x) G{x + a, x + a). Similarly the bubble instability-the fact that 
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FIG. 5: The self-consistent bounce action for both models is plotted as a function of the bounce 
radius as the bounce is dilated and contracted. The curve fits are quadratic functions and require 
only the location of the peak and the second derivative (curvature) as parameters. 

it is a saddle point of the effective action-can be probed via a simultaneous radial dilation 
or contraction of both (j) and G, see Fig. |S] 

Care must be taken when tracing over the collective coordinates because the contribu- 
tions from these degrees of freedom carry dimensions and potentially lead to new (infrared) 
divergences. The regularization of such behavior should not introduce new ad hoc scales in 
the calculation but rather should be consistent with the regularization and renormalization 
implemented in Sec. IIVI In practice, the dimensionality introduced by the collective coor- 
dinates is compensated for in the transition amplitude by an infrared regulator /xq, which 
is tied to the choice of renormalization scale made above in dealing with the logarithmic 

n 

divergences in the effective action 31]. The physical (measured) mass of quasiparticles is a 
function of this scale M = M{fio). Thus by fixing the mass in the false vacuum, we have 
implicitly adopted the scale choice /xq = ^. 

^ Another possible scale is given by the inverse object radius, i.e. = ^ E^- As can be shown 



We now evaluate the (consistently regularized) contribution from the collective coordi- 
nates, beginning with the translational zero modes. The effective action is invariant under 
a transformation which shifts the center of the bounce. As usual, the path integral obtains 
a factor of spacetime volume up to a Jacobian factor Jq, which in D dimensions is given by 



D[(l){x, z)] = Jod^z^ with Jq 
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— / d^'x 



D / o<Pb 



2tt \ dz 



1 D/2 



(69) 



where the factor of 271 is conventional to make the integral commensurate with a the Gaussian 
case. 

The integral can be expressed in terms of the effective action. To be explicit, we wish to 
consider the effective action now as a functional of alone. This amounts to the familiar 
transformation of the 2PI action rcfr[0, G] into the IPI effective action rgg^[0], wherein G is 
thought of as a functional of determined by Eq. (fTTjl jl^. Dividing this action into kinetic 
and effective potential terms, the equation is simply 

- V V + ^ = 0. (70) 

and we have essentially the same situation as in the 1-loop analysis but for the substitution 
of an effective potential. The argument proceeds that the effective action is conserved along 
a trajectory in x: dxT^s = and its value is fixed by our choice of renormalization in the 
false vacuum reff[0+] = 0. As a consequence we can write 



1 f d(j)Y T. n T [reffl^'^^ ^^^^ 



2 \dx J 
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The total contribution from the zero modes, together with factors of the regulator /iq, is 

Jof^o [ d"^- (72) 



Next we consider what happens upon variation of the radius corresponding to the bounce 
solution in and G. It is convenient to shift the radial coordinate by the bounce radius, 
i.e. r = i?o + ^, where r = Rq is the radius of the bounce. Expanding the effective action 
around ^ = then gives 

Teff ^ refr(i?o) + iK^{Ro) + \eKs{R^) + ■■■ (73) 



explicitly in the thin wall approximation, this latter choice ties the regulator scale to the symmetry 
breaking parameter e. 



where primes denote radial derivatives. The first derivative of the effective action at the 
bounce is zero as the bounce is a solution of the equations of motion. The second derivative 
however is negative, as shown in Fig. El Just as is the case for the semiclassical solution, 
the self-consistent bounce is a saddle point of the effective action and a maximum in the 
radial "direction" (by which we mean direction in function space). The contribution of this 
dilatational degree of freedom to Z is computed in direct analogy to that of the negative 
mode in the 1-loop analysis, i.e. in steepest descent by continuation to imaginary ^ jj]. 
The difference is that we must now extract F^^'g by hand, so to speak, rather than from the 
fluctuation spectrum. Performing the integration in the radial coordinate requires again the 
introduction of a Jacobian 
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(74) 



The change to radial coordinates has no effect, and J_ is identical to Jq ^oi one dimension. 
The -D + 1 collective coordinates give a total multiplicative contribution to Z of 
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(75) 



or, performing the Gaussian integration. 



. D+1 

27r 



27r 

^ eff 



(76) 



where the factor of 1/2 results from the analytic continuation to one branch of imaginary ^ 
since F^g is negative. 

Finally the rate of spontaneous decay of the false vacuum is obtained by Eq. (jH)), 



T/V = fio 
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(77) 



where all factors of Feg and its derivatives are evaluated at the self-consistent bounce. The 
dimensions are correct as the factor containing derivatives of the effective action cancels one 
mass dimension. We present numerical results for the nucleation rate computed according 
to our self-consistent method described above and computed at 1-loop order in Table |nj 
The models are as in Eqs. ()18p - (fT^ with parameters given by the numerical values /i^ = 



4.34, A = 0.14, e = 0.67 for the 6'^ model and A 



mu 



^ = 1, e = 0.12 for the (b'^ model. 



(f)^ model 1.19 x 10""^ 8.59 x 10"^ 
(f>^ model 1.23 x 10"^ 1.60 x 10"^ 

TABLE II: Numerical results for the self-consistent nucleation rate and the 1-loop rate in the two 
scalar models. 

VII. DISCUSSION 

We see that the self-consistent calculation results in a much smaller nucleation rate than 
the usual 1-loop estimate, from two orders of magnitude in the 0*^ model up to six orders of 
magnitude for the 0^ model. This is a robust prediction even if the nonexponential prefactor 
in Eq. ()77|) is taken with a grain of salt. The prefactor is of the same order of magnitude 
in both the self-consistent and 1-loop approximations, hence it is the exponentiation of 
the effective action which mostly determines the magnitude of the rate. We should clarify 
that the "dynamical" prefactor, i.e. the ratio of fluctuation determinants in Eq. (^, is not 
generically of order unity. It is only the case that when the (appropriately dimensionless) 
determinant contribution is absorbed into the effective action in the exponent that it may 
be argued that the remaining prefactor is roughly determined by some mass scale typical of 
the problem as . 

The larger effective action at the self-consistent level results principally from the the shift 
in the spectrum (relative to 1-loop) by the fluctuation self-repulsion. The same qualitative 
change is observed for self-consistent topological defects, which are heavier than at 1-loop 
[lfil |. and we expect it to be generically true for any self-consistently dressed quantum field 
configuration. 

In summary, we have proposed an alternative nucleation rate computation based on a 
self-consistent determination of the bubble shape and fluctuation spectrum. We have traced 
a path analogous to the standard work of Coleman ^3] by going to imaginary time and 
solving instead a set of equations for the mean-field bounce and coincident two-point function 
obtained from a 2PI effective action. As translational invariance and bubble instability are 
not manifest at the level of the linearized fluctuation spectrum, we have used collective 
coordinates to account for the contribution of these degrees of freedom to the transition 
amplitude. 




FIG. 6: A stroboscopic view of the bubble growth for the self-consistent bubble of the (f)^ model. 
The bubbles evolve by expansion of the true vacuum and Lorentz contraction of the walls, without 
giving rise to radiation, as required of the self-consistent bounce. 



Finally we note that as our Euclidean spacetime solution, upon analytically continuation 
to Minkowski spacetime should provide a solution to the real time bubble dynamics for all 



221. This is the ultimate 



time in accordance with Lorentz invariance at zero temperature 
test of whether or not what we claim above to be the self-consistent bounce (at this level of 
approximation) is indeed the self-consistent bounce. 

Analytic continuation of the one- and two-point functions 0(p) and G{p) from the last 



section is a straightforward matter of taking p —>■ \Jx^ — t^, or at initial time t = 0, p — > 
This is not however enough information to specify the real time problem. Time evolution 
of the coupled system of equations requires full knowledge of the spectral decomposition of 
G[x^x'\ In order to initialize these, we need to rebuild in real time the fluctuation modes 
that correspond to the (initially static) 0(x) and G{x,x). In practice we use the methods 
developed in Ref. ^| to easily determine the mode functions that correspond to </)(x) and 
G(x, x) found in the previous sections. We then use them as initial conditions for the real 
time field evolution. 




At zero temperature, we should therefore observe a Lorentz invariant bubble growth 
without the initial transients present in the naive initializations of Ref. 22]. As figures IHl 
andEl corresponding to the 0^ and (j)'^ models, indicate the shape of the bubble is essentially 
unaltered during the evolution, except for Lorentz contraction. 
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